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<N . ABSTRACT 

a\ 

^ ■ The afterglows of gamma-ray bursts (GRBs) have commonly been assumed 

to be due to shocks sweeping up the circum-stellar medium. However, most 
GRBs have been found in dense star-forming regions where a significant fraction 
of the prompt X-ray emission can be scattered by dust grains. Here we revisit 

1 the behavior of dust scattering of X-rays in GRBs. We find that the features 

of some X-ray afterglows from minutes to days after the gamma-ray triggers 

' are consistent with the scattering of prompt X-ray emission from GRBs off host 

dust grains. This implies that some of the observed X-ray afterglows (especially 
those without sharp rising and decaying flares) could be understood with a dust- 

■ scattering-driven emission model. 

Subject headings: dust, extinction — gamma rays: bursts — interstellar medium 
— X-rays: general 

V , 

1. INTRODUCTION 

a ; 

Currently, the popular model for gamma-ray bursts (GRBs) and their afterglows is the 
fireball-shock model (for recent reviews, see Zhang & Meszaros 2004; Piran 2005; Meszaros 
2006). In this model, the short-term prompt emission of GRBs is ascribed to internal shocks 
in the ejecta, and the long-term afterglow at lower energy is ascribed to external shocks 
sweeping up th e circum-stellar medi um. The rapid localization of the Swift Gamma- Ray 



Burst Explorer ( iGehrels et al. 1120041 ) has led to a recent breakthrough in the detections of 
early afterglows. This has given rise to theoretical studies of the early afterglow emission, 
especially in X-ray energy, thanks to the X-Ray Telescope (XRT) onboard Swift. 



Recently, an X-ray halo around the short GRB 050724 was detected by XRT, with both a 
radial temporal evolution and intensity distribut ion that are consistent with the properties of 
Galactic dust-scattering (jVaughan et al. II2006I ). Previously, GRB 031203 (IVaughan et al. 
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2004 ) and GRB 050713 A (ITiengo fc Mereghetti 1120061 ) had also been found to have similar 
dust-scattered X-ray halos, as observed by XMM-Newton. In addition, a careful inspection 
of the XMM-Newton data shows the presence of some diffuse emission in the field of GRB 
050730, which could also be produced by dust-scattering, although this needs to be verified 
( jTiengo fc Mereghetti 1 120061 ) . Note that host galaxies have been found for some GRBs 
localized so far. Most of these galaxies show signs of active star formation, imply ing the 
presence of GRB progenitors forming out of dense gaseous clouds ( iPaczyriski Ill998l ). Thus, 
dust-scattering off of dust grains may be common in GRB phenomena, which could have 
been playing an important role in the observed X-ray afterglows. 

A diffuse X-ray halo is predicted to appear around an X-ray po int source when the inter- 



stella r dust grains scatter some of the X-rays, typically by 1' to 1° ( lOverbeck Ill965l ; iMartin 



197oh . This time-depend ent information about scattering in GRBs was previously considered 



by lDermer et al. I (119911 ). who assumed the existence of binary companions or accretion disks 



in GRB systems; a Compton echo of reflected X-ray and gamma-ray emission with a time 



emission were also discussed in detail by 


Miralda-Escude 


f 


1999), 


Meszaros & Gruzinov 


(2000h. 


Madau et al. 


f 


2oooh. 


Esin & Blandford 


f 


2000h. f 


Sazonov & Sunvaev ( 


20031) and 


Ramirez-Ruiz & Madau 


( 


2004 


) in a variety of emission geometries and ambient gas distri- 



butions around GRBs. 

In this paper, the scattering of X-rays by dust grains in GRBs is revisited. We estimate 
the emerging flux during such an X-ray echo event, which is expected to be dominant in the 
X-ray afterglow. We find that an initial pulse of X-rays from a normal GRB scattering off 
dust grains in a host galaxy can give rise to a long-term "afterglow" with almost the same 
amount of energy as expected in GRB X-ray afterglows. However, the angular size of this 
echo emission is too small to be possibly resolved on Earth. Therefore, only the temporal 
features of the total scattered flux are considered here. We find that some of these features 
are consistent with observations of GRB X-ray afterglows. We discuss the prominent flux of 
X-ray echo emission in § 2 and the temporal behavior in § 3. We suggest in § 4 that dust 
scattering of prompt X-ray emission off host dust grains may be an alternative explanation 
for some of the GRB X-ray afterglows. Our conclusions are summarized in § 5. 



2. TOTAL FLUX OF AN X-RAY ECHO 

The quantity that we first consider is the amount of energy. This should be evaluated for 
the scattering effect in the context of GRBs at the first step to make sure that scattering off 
dust grains can viably produce a detectable X-ray afterglow. This was previously discussed 
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by iMiralda-Escude I (119991 ) . We follow his result here. 



Considering a variable X-ray source at angular diameter distance D s , and an intervening 
dust layer (e.g. a galaxy or a cloud) at redshift Zd and with angular diameter distance D d , 
the time delay t d of photons received at an angle 9 from the source scattered by the dust 
layer is given by 

. (l + z d )D d D s 9 2 
2cD ds 

where D ds is the angular diameter distance from the dust layer to the source. See Figure 1 
for the geometry of the scattering. 

For simplicity, the variable source is assumed to emit most of its energy in a narrow 
energy band (e.g. 2-10 keV) with a given ffuence So, in a time much shorter than t d . Thus 
we can approximate the source emission as a pulse of monochromatic light with wavelength 
A. If an experiential size distribution of grains in the intervening dust layer is assumed, i.e., 
r(a) oc a -0,5 , where r(a) is the scattering optical depth of dust grains with radius a, the 
total flux of the X-ray echo due to the small-angle scattering can be estimated by 

F h ~ 1.96 x 10- 7 ° ' p in22 b , (2) 

D ds t d X 10 22 cm 2 0.02 

where a p ~ 1.5XD ds / '(2ttD s 9) is approximately satisfied, N b is the mean gas column density 
over the layer, and Z t is the total metallicity in dust grains. Here the differential cross 
section is treated as an approximation of the Mie solution, which is given in the form of the 
"Rayleigh-Gans" approximation in § 3. 

Using the relation between t d and 9 (eq. [1]) in equation (2) yields 

Fh[td) — 48.5 ^ — j A 3/ 2 10 22 cm -2^> ( 3 ) 

where the echo flux is a function of the time delay t d . We should point out that (1) this 
equation is a good estimate of the flux of X-ray echoes, from which we can get a preliminary 
idea of how much the flux will be at a given t d and whether it will be detectable at the 
time of interest. We discuss this issue later in this section. (2) Considering the typical 
scattering angle implied by the value of a p , equation (3) is only valid before or around 
t d ~ 10 10 s(9A 2 /32tt 2 c)(1 + z d ){DdsD d /D s )(a p /0.1 fim)- 2 , where c is the speed of light. (3) 
This equation implies that a pulse of light propagating through a dusty region will produce 
a long-term delayed emission, which is like an "afterglow" of this pulse (see Fig. 1). Here 
we call this a dust-scattering-driven afterglow. We work out its detailed light curve in § 3. 

Equation (3) suggests that the flux of echoes is substantially determined by the given 



-4- 



geometry of dust scattering (i.e., D s , Dd, and Dds) at a given time. In general, three cases 
should be considered. 

Case (1): the dust layer is in the form of an intervening galaxy, which is also at a 
cosmological distance. Thus, 

x 3/2 / +. \ -V4 

V 



F h 



1.27 x 10" 15 
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10 22 cm" 2 0.02 



10 _6 ergs cm -2 
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6keV, 



td_ 
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(4) 



where e is the mean energy of the X-ray photons, -D<f/(1 + Zd) — 10 3 Mpc, and D s /Dd s — 2 
are assumed. This X-ray echo can be mu ch brighter than any intrinsic emission from the 
intervening galaxy (IMiralda-Escude Ill999l ). However it is belo w the XRT detectio n limit 
(~2x 10~ 14 ergs cm -2 s" 1 for a 10 4 s integration time) of Swift (IGehrels et al. II2004I ). 

Case (2): the dust scattering happens at low latitudes or in dense clouds of our Galaxy, 
which have been verified by observations. We can get a similar expression for the total flux 
of the echo 
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wher e Zd 



0,D S ~ 



Ddsi Da 



100 pc, and N b 
19901 1. The typical halo radius is about 6 



10 21 cm 2 are assumed (IDickey fc Lockman 
1.5'(t (i /10 3 s) 1 / 2 (D d /100pc)- 1 / 2 . This echo 
feature is quite detectable and consistent with observations (Vaughan et al. 2004, 2006). 
Due to the large halo radius , this case can be easily distinguished by observational analysis 
jTiengo fc Mereghetti Il2006h . 



Case (3): if the dust scattering happens in the host galaxy of a GRB at redshift Zd 



we have 



2.25 x 10~ 9 

/ Dds 
VIOOpc 
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(6) 



where D s ~ Dd, Dd s — 100 pc, and N ~ 10 22 cm -2 are assumed. This case is similar to case 
(2), and thus the echo emission is significantly detectable, since both cases have a similar 
scattering geometry, except that the halo in case (3) has a very small angular size and cannot 
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be resolved on Earth. Obviously such a bright flux is quite comparable to those of normal 
X-ray afterglows (e.g. Costa 1999). We need to give a further consideration of this echo 
component, which may have been observed but not realized so far. In general, only the total 
flux of the echo emission versus time needs to be considered, due to its small angular size. 



LIGHT CURVE OF AN X-RAY ECHO 



Here, we revisit the temporal behavior of an X-ray echo event with a delayed time 
of minutes to days (e.g. Meszaros & Gruzinov 2000; Sazonov & Sunyaev 2003). For this 
purpose, we need to know the differential cross section for the small-angle scattering of X-rays 
off dust grains. For typical spherical grains of radius a, in the limit of 



IkeV 



1/rni 



< 1 



(7) 



the exact Mie solution recovers the Rayleigh-Gans approximation 



Xr 



where a? is the total cross section, A is the wavelength of the X-ray photons, x = (2na/\)a 
is the scaled angular coordinate, a is the scattering angle, and h(x) = (sinx)/a: 2 — (cosx)/x is 



the fi rst-order spherical Bessel function (lOverbeck II 1 965l ; I Alco ck fc Hatchett Hl978l ; lSmith fc Dwek 
19981 ). For simplicity, our treatment takes the dust grains as Rayleigh-Gans particles (e.g. 
Kriigel 2003) and is independent of the chemical composition or shape of the grains, as long 
as equation (7) holds. 



3.1. Analytical Treatment in a Simple Case 

For a variable X-ray source (e.g. the prompt X-ray counterpart of a GRB) at an angular 
diameter distance D s , with an unabsorbed flux as a function of time, F u (t), scattered by an 
intervening dust layer (which could be assumed to be in the host galaxy of a GRB with an 
angular diameter distance Dd s from the GRB) at an angular diameter distance Dd with a 
scattering optical depth t(6, 0), the observed intensity / of the X-ray echo is calculated by 

I(e, * t) = F.( t - ti )r { 0, *)-%-^- = 4 " D - K(t -*' )T f' . 0) 

where td is given by equation (1), a = (D s / D ds )9, dflsc — adadcj), and dQ = 9d9d(f> are 
used in the small-angle limit, A is the mean wavelength of X-ray photons in the dust frame, 
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r(9, <ft) is assumed to be small, and the light that multiply scatters into the line of sight is 
ignored. 

At first, for simplicity, we assume that r(6, <fi) = r does not vary with 9 and 4> based 
on the fact that the dust layer may have a small thickness (Vaughan et al. 2004, 2006). For 
a bursting source, most of its energy is assumed to be emitted in a narrow band on a time 
much shorter than t d , and the flux is then approximated as F u (t) = Sod(t), where 5(t) is the 
Dirac delta function and the source's trigger time is taken as the time zero point. It follows 
from equation (9) that the total flux of the echo is analytically expressed as 



F h (t) = / I(9,<P;t)cos9dtt 



^S aD s T aJ l{x[6{t)}}6{t) 
XD ds x[9(t))t 



(10) 



where the function 9(t) is defined as 9{t) = [2ctDd s /(l + Zd)DdD s } 1 ^ 2 , and x{9) is defined as 
x(9) = 2naD s 9/(XDds)- Some formulae are given in the Appendix for translating the Dirac 
delta function of t into a function of 9. 

This temporal behavior is illustrated in Figure 2 with the dot-dashed line. The dust 
layer is assumed to be in the host galaxy, and we have Dd = D s , Zd — 1, and Dd s = 100 pc in 
Figure 2. Studies of interstellar extinction indica te that most of the grains have a si ze near 
a ~ 0.1 /im and a distribution in a wide range (IMathis et al. 1 119771 ; iDraine 1 120031 ). Here 

we choose a = 0.1 /xm for a simple illustration. Th e X-rays in the range 2 6keV can 

be scattered most efficiently (IMiralda-Escude Ill999l ). and thus the average photon energy 
e = 6keV is used. The total X-ray fluence Sq ~ 10~ 6 ergscm -2 in the energy band from 
about 2 to 10 keV is given. The mean scattering optical depth To = 0.1 is assumed in Figure 
2, although r is expected to be dependent on the photon energy (e. g. Meszaros & Gruzinov 
2000; we account for this below). 

Note that in the X-ray scattering scenario, the size of the interstellar dust grains is 
generally larger than the wavelength of the X-rays. As very big particles block practically all 
light that falls onto them, most of the scattered photons substantially originate from diffrac- 
tion at their edges. The diffraction of X-rays through a dust layer causes an interferometric 
pattern in the differential cross section as a function of the scattering angle, as given in 
equation (8), only if equation (7) holds (e.g. Kriigel 2003). This translates into the pattern 
of total flux versus time, as shown in Figure 2. The significant semi-periodic interferometric 
pattern indicated in these light curves reminds one of the X-ray light curve o bserved in some 
long GRBs, e.g., GRB 050904 flWatson et al. Ibood Icusumano et al. Il2006h . Of course, this 
pattern may be smoothed out in some other realistic situation, when there is a range of dust 
grain sizes and the total flux is detected in a finite X-ray band (e.g. 0.3 — 10 keV for XRT 
). This is treated in some detail in the following subsection. 
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The shallow decay phase shown in the light curve, which is common in the small- 
angle scattering scenario predicted by equation (3), can actually be attributed to the first 
maximum in the differential cross section peaking at a very s mall scattering angle, wit h 



1978). 



1978|) 



the scaled notation x = (27ra/A)a ~ 2naD s 8/(XD ds ) ~ 1.5 flAlcock fc Hatchett 
The differential cross section decreases dramatically at x > 3.0 (lAlcock fc Hatchett 
which is translated into the decaying total flux at t > (1 + Zd)DdDd s \ 2 / (SiTcD s a 2 ) ~ 
5 x 10 3 s(e/6keV) _2 (a/0.1/xm) _2 (.D (is /100pc). After that, a fast power-law decay of the 
maxima (roughly oc t~ 2 ) emerges. Regardless, the small-angle approximation used here 
(i.e. sin a ~ a) is not violated throughout the whole light curve, since the scattering angle 
a ~ (ctd/ Dds) 1 / 2 is no more than several arcminutes at a delayed time — 10 5 s, with 
D ds ^ 100 pc. 

Furthermore, we also consider the scattering of a constant X-ray beam emitted by the 
source within a time range T, i.e., the flux is approximated as F u (t;t < T) = Fq. The 
received flux can be calculated by integrating equation (9) over 8 and (ft. The results with 
different T are also shown in Figure 2. The on-going flux Fq = 10~ 6 ergscm~ 2 s" 1 is given, 
and the other parameters are not changed. The light curves due to scattering of a beam are 
similar to those of a pulse, except that an early rise is expected in the former case. 



3.2. Numerical Treatment for a Practical Case 



Above we assume a constant scattering optical depth tq, a grain size a, and a homochro- 
mous initial fluence So. In some cases , r may vary as the X-ray energy e changes, and the 
grain size may have a certain distribution. For the Rayleigh-Gans approximation, the total 
scattering cross section <r scal 



if. a) oc e 2 a 4 is also inferre d (e.g. Mauche fc Go renstein 1986), 



and it is suggested by observations that n(a) oc a q (IMathis et al. II1977I ). Here, in our 



treatment with only small scattering optical depth (r oc na sca ), we assume that 



r(e, a) = A 



IkeV 



0.1/iHl 



4-g 



for 



0.2 < e < lOkeV ' 



where A = (5 — g)rk e v/[o(a/0.1//m) q }\^+ is a constant in units of cm , Tk e v is the scattering 
optical depth at 1 keV, q ^ 3.5 — 4.5, a~ — 0.005 — 0.025 Mm, and a+ — 0-25 — 0.5 ixm ar e 



inferred from observations (IMathis et al. 1119771 ; iMauche &: Gorenstein Ill986l ; iDraine II2003I ). 
Here s ~ 2 is inferred from a single observation ( Mitsuda et al. 199ol ) and then adopted by 
Meszaros & Gruzinov (2000). 



An initial source spectrum should also be taken into account. In the soft X-ray band, 
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this is given by the Band spectrum (IBand et al. 1119931 ) 

S(e) = B ( 



.lOOkeV, 





r (* + l)ej 


1 exp 





(12) 



where B is a parameter in units of ergs cm -2 keV -1 , and 5 ~ and E p ~ 200 keV are 
suggested by Preece et al. (2000). Thus, we can calculate the measured flux, e.g., in the 
0.3 — 10 keV band by XRT. Consequently, 



Fk(t) 




47rABaD s T(e } a)jf{x[e(t)]}§(t) 



\D ds x[6(t)]t VlOOkeV 







) exp 





dade (13) 



Here the unabsorbed initial emission is still taken to be a pulse of light, approximately as 
in equation (10). In the case of a beam of light, the difference in the echo light curve is the 
short rising time at the beginning time (e.g., see Fig. 2). 

Equation (13) can be evaluated numerically with different parameters considered. We 
find that the temporal behavior of significantly depends on three observable quantities 
of dust grains. The first quantity is the position of dust D ds . This is because t d depends 
prominently on the position of the dust, which is inferred from equation (1). The second 
quantity is the maximal size of a dust grain, a + , which is the large-size cutoff for the size 
distribution. The third quantity is the index s in the dependence of dust-scattering optical 
depth on X-ray energy. 

One prominent characteristic of this temporal behavior is a distinguishable broken power 
law, with a break time dependent on the parameters discussed above. As shown in Figure 3, 
a shallow decay before the break and a slope ~ —2.0 after the break are clearly present in the 
light curves. The shallow decay before the break is expected from our previous discussion 
and roughly consistent with our estimates of small-angle scattering, shown by equation (3), 
while the steep decay after the break can be attributed to the quickly decreasing cross section 
of larger angle scattering at later times. This feature of the flux (oc t~ 2 ) is also implied in 
the previous discussion, e.g., as shown in Figure 2, which is roughly consistent with the 
decreasing maxima. 

Figure 4 plots the spectral evolution during dust scattering. In general, the spectra 
soften as the flux decreases. The softening of the spectra can be attributed physically to 
the diffraction effect (treated as scattering here) which is described in the differential cross 
section versus scattering angle (eq. [8]). Softer X-rays tend to be scattered at a larger 
angle a, with high-order maxima in the differential cross section, and are thus received at a 
larger angle 6, which leads to a longer arrival time due to a longer light distance (eq. [1]). 
Obviously, one can tell from Figure 4 that the visible softening emerges at later time when 
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the steep decay (oc t 2 ) in the light curves shows up. This corresponds to a delayed time of 
t d ~ (l + z d )D d D ds X 2 /(32tt cD s a 2 ) ~ 3 x 10S(e/lkeV)- 2 (a/0.1/mi)- 2 (Ai s /100pc), where h is 
the Planck constant, and A, e, and a are the equivalent average photon wavelength, energy, 
and dust radius, respectively 



IMPLICATIONS FOR X-RAY AFTERGLOWS 



Above, we revisit the X-ray scattering off dust grains at GRB stages. Now we apply 
these results to observational data and suggest that some of the X-ray afterglows detected so 
far may be alternatively explained as the emission from X-ray echoes, i.e., dust-scattering- 
driven afterglows. 

1. A shallow decay followed by a "normal" decay and a further steepening is suggested 
by Zhang et al. (2006) to be characteristic of almost all the Swift GRBs. (1) In gen- 
eral, to account for shallow decay, a continuous activity of the GRB progenitor is expected 
( IDai &: Lu I Il998t IZhang. fc Meszaros I l200ll; iDai I I2004T) or a pow e r law distribution of th e 
Lorentz factors in the ejecta is assumed faees &: Meszaros" 1998 ; Sari fc Meszaros 2000 ). 
Alternatively, we propose that this feature can be attributed to the X-ray echo emission at 
early times. (2) To account f or the steep d e cay after the "normal" decay , the relativistic 
jet effect has been suggested ( IRhoads Ill999l ; ISari. Piran. &: Halpern Ill999l ). Alternatively, 
this feature can be explained as the X-ray echo emission at late times. To summarize, the 
features with a shallow decay followed by a "normal" decay and a further steepening are 
consistent with the X-ray echo emission presented above, e. g., as shown in Figure 3. 

Here, we apply equation (13) to two recently detected GRBs, 060813 and 060814, in 
Figure 5, where 5^0 and E p ~ 200 keV are assumed. The consequent parameters are 
a + ~ 0.5 /mi and D ds ~ 10 pc for GRB 060813 and a + ~ 0.25 /mi and D ds ~ 30 pc for 
GRB 060814. The other parameters, a_ ~ 0.025 /mi, q ~ 4, and s ~ 2, are the same for the 
two GRBs. Obviously, a shallow decay is common in early X-ray afterglows. This favors the 
dust-scattering scenario, which predicts a shallow decay at an early time, when the scattering 
angle is smaller (see also eq. [3]). In any case, the detailed light curve depends on several 
parameters that we mentioned in §3.2, and thus the early temporal index varies in a wide 
range, e.g., ~ [0, —1.0], and then steepens into ~ [—1-0, —1.5] at a moderate time (see also 
Fig. 3). Of course, as the scattering angle gets larger at later times, a further steepening of 
the decay (oc t~ 2 ) is also predicted in our c alculations. Never t heless, the steep decay (oc t~ 2 ) 



2007; Romano et al. 



2006bf ). In the 



is not observed in most X-ray afterglows flSato et al. 
latter case, it may be due to a larger D ds , which means that an echo event with a small 
scattering angle will take place with a longer duration, so that the shallow decay component 
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will last a longer time, thereby preventing the emergence of a steep decay. 



2. Most early X-ray light curves are fou nd to decline rapidly i n the first few min- 



utes, with a power- law index of ~ 3 or greater (jTagliaferri et al. 1120051 ; iNousek et al. 112006 



Zhang et al. 1 120061 ) . In gen eral, this feature can be tak en as a GRB tail emission arising 
from high angular latitudes (IKumar fc Panaitescu 1 120001 ) . Then a break to a shallower de- 
cay com ponent, as men t ioned above, commonly shows up, which is defined as an X-ray 
hump by lO'Brien et al. I (120061 ). Early XRT data reveal little evidence of spectral evolution 
across this temporal b reak, but where evo lution is seen, the spectrum tends to get harder 



( jO'Brien et al. 



20061 ; INousek et al. II2006I ). Provided that the "hump" component is ex- 



plained as the X-ray echo emission, this spectral feature is basically consistent with our 
expectation. As shown in Figure 4, the early-time spectrum of echo emission in soft X-rays 
would hardly change from the initial prompt spectrum. Given an intrinsic evolution of the 
prompt spectrum, the delayed echo spectrum can be harder than the contemporaneous tail 
emission in the observer's frame. 

Our model predicts that a visible spectral evolution in XRT is expected when the steep 
decay phase (oc t~ 2 ) emerges. It is worth mentioning that the Rayleigh-Gans approximation 
is adopted here in calculating the differential cross section. This approximation works pretty 
well for normal interstellar dust an d energies at or above 2 keV. However, it overestimates 
the echo emission at lower energies (ISmith fc Dwek Ill998l ). where absorption of soft X-rays 
is important. This absorption effect will weaken the low-energy component in the spectra. 
Thus, the spectral shapes, as shown in Figure 4, have been idealized in our treatment. 
Regardless, a softening spectrum is expected when the steep decay phase (oc t~ 2 ) emerges. 
This feature needs to be verified by time-resolved spectral analysis, in case the steep decay 
phase does emerge. 

3. Fluctuations or flares are observed in some GRBs (e.g. GRB 050904: Watson et al. 
2006, Cusumano et al. 2006; XRF 050406: Romano et al. 2006a; GRB 060713A: Guetta et 
al. 2006). Generally these flares are thought to be caused by late internal shocks, similar 
to those that produce the prompt emission (Zhang et al. 2006; Burrows et al. 2005b; Fan 
fc Wei 2005; Wu et al. 2005; Perna et al. 2005; Liang et al. 2006), or late external shocks 
( IPiro et al. II2005I ). Both the late inter nal shock model and late external shock model require 
late-time activities of central engines (jWu et al. 1120051 ). 



Alternatively, we have three reasons suggesting that these flares, or some of them, might 
be due to dust-scattering: 

1. Enormous flares are also observed in some GRBs (e.g. GRB 050502B: Burrows et 
al. 2005a; Falcone et al. 2006), which requires a more complicated theoretical explanation. 
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However, these flares usually begin with an excess absorbing column that softens as the flare 
progr esses. In addition, the afterglow intensity a nd slope are similar before and after the 
flare (IBurrows et al. I l2005al ; iFalcone et al. 1 120061 ) . These features seem to favor the dust- 
scattering scenario, which does not disturb anything when it turns off. The softening feature 
emerges when it turns on. 

2. Semi-periodic fluctuations are revealed in some GRBs (e.g. GRB 050904: Watson 
et al. 2006, Cusumano et al. 2006; GRB 050730: Burrows et al. 2005b), which reminds 
one of the interferometric pattern expected in the dust-scattering scenario under certain 
assumptions. Although this feature may always be smoothed out in practical cases (see § 
3), one can see a hint of the physics this may have unearthed. 

3. Otherwise, if the dust layer is broken up or distorted by the GRB progenitor or is 



fluffy, as in some cases with many holes and voids (I Woo et al. 1 1 19941 : iMathis et al. 1 11995 



Predehl fc Klose II 19961 ). flares should also be expected in the echo light curves. For example, 
a big jump in r{6) around 6q causes a rapid rise and a rapid decay around t d — (1 + 
z d )D d D s 6l/(2cD ds ) ~ 10 3 s((l + z d )/2)(£> dfl /100pc)(ao/60'') 2 due to the time delay of the 
scattered flux introduced in equation (9), where a is the angular scale from the GRB source 
of the fluctuation of the dust intensity in the host galaxy. In any case, we may be able 
to discover whether these flares are due to dust scattering, after improving our model with 
specific dust details. 



5. CONCLUSIONS 

In this paper, we revisit the X-ray dust scattering in GRBs. First, we give an estimate 
of the emerging flux during such an X-ray echo event, which is expected to be dominant in 
X-ray afterglows. We find that an initial pulse of X-rays from a normal GRB scattering off 
dust grains in the host galaxy can be followed by a long-term "afterglow" with almost the 
same amount of energy as expected in the GRB X-ray afterglows. Second, we investigate the 
behaviors of such an event, especially its light curves. We find that some of these features 
are also consistent with observations of GRB X-ray afterglows. We suggest that some of the 
X-ray afterglows from GRBs (especially those without sharp rising and decaying flares) can 
be understood in the dust-scattering-driven afterglow model. The scattering of the prompt 
X-ray emission from GRBs off the host dusty regions can be an alternative explanation for 
most of the features observed recently in X-ray afterglows by Swift XRT. 

Several properties of dust grains (e. g., D ds , a+, and s) are supposed to be relevant 
for temporal behaviors of X-ray echoes. The most deterministic one is the position of dust 
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(i.e. Dds)-, because it determines the time delay td given in equation (1). However, until 
now, we have not known this quantity very well. Reichart (2001) assumed that there is 
a preburst, dense environment due to the strong winds of GRB progenitors with an inner 
radius of about several parsecs, while Madau et al. (2000) and Moran & Reichart (2005) 
suggested that this radius may be about 0.001-0.01 pc. Meszaros & Gruzinov (2000) and 
Sazonov & Sunyaev (2003) assumed a GRB origin at the center of a uniform dusty region 
with a radius of about 10-100 pc. In addition, the situation could be more complicated, with 
the evolution of the dust grain population considered. Waxman & Draine (2000) suggested 
that dust grains will be sublimated by the optical/UV flash of GRBs out to a distance of 
about several parsecs (see also Perna & Loeb 1998; Lazzati et al. 2001; Perna & Lazzati 
2002; Heng et al. 2007). Recently, Campana et al. (2007) suggested a distance of several 
parsecs by analyzing the evolution of the soft X-ray absorbing column around GRB 050904. 
Here, we find that a D^s ~ of ~ tens of parsecs is consistent with Swift GRBs, based on our 
dust-scattering-driven afterglow model. 

It should be noted that X-ray echo emission is proposed here to be due to dust scattering 
taking place at a distance of ~tens of parsecs from the GRB source. However, in the standard 
external-shock model, relativistic shocks generally take place up to a distance of ~ 10 17 cm. 
So technically this echo scenario does not rule out the existence of emission from external 
shocks in both X-ray and optical/NIR bands. In addition, there is indeed some evidence for 
chromatic light-c urve breaks, which may require that the X-ray and optical emissi on have 
different origins (IFan fc Piran 1120061 ; iPanaitescu et al. 1120061 ; iRomano et al.l l2006bl ) . Thus, 
there will be at least two types of X-ray afterglows if the external-shock model and the 
dust-scattering model are both valid. 
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A. APPENDIX - FORMULAE FOR ANALYTICAL TREATMENT 

The Dirac delta function has the following properties: 

6(-x) = 5{x) (Al) 

8(ax) = |a|- x 5(x) (A2) 
S( x 2 - a 2 ) = (2\x\)- 1 [5(x + a) + 5(x - a)] (A3) 

In § 3.1, F u {t — t d ) = S 5(t — t d ) is assumed for a GRB pulse. Here, the delta function 
of t can be translated into a function of 9: 

2cD ds rr/l2 2ctD ds . 



[l + z d )D d D. 1 (l + z d )D d D 
cD ds 



1 + z d )D d D s 9 

cDds 
1 + z d )D d D s 9 



{6[9 + 9(t)) + 5[9-9(t))} 
8[9 - 9(t)} (A4) 



where the function 9(t) is defined as 9{t) = [2ctD ds /(l + z^D^JJ ^l 2 , and 9 is always positive 
in our treatment. 
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Fig. 1. — Scheme of X-ray small-angle scattering. The scattering angle is magnified to be 
illustrative. In fact, the dust layer must be very close to the line of sight. 
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Fig. 2. — Light curves of dust-scattered X-ray echoes. T is the timescale of the X-ray 
beam defined in the text; the pure pulse is T = 0. Parameters So = 10 _6 ergscm -2 , Fq = 
10~ 6 ergscm~ 2 s" 1 , D s = D d , Zd = 1, D ds = lOOpc, a = 0.1/xm, e = 6keV, and r = 0.1 are 
assumed. 
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Fig. 3. — Light curves of dust-scattered X-ray echoes with detailed dust properties con- 
sidered. The solid line show the model with normal parameters a_ = 0.025yum, a + = 
0.25/im, q = 3.5,Dd. s = lOOpc (we assume that the dust layer is in the host galaxy), and 
s = 2.0. Otherwise, we let only one parameter change in each line. ^From top to bottom, 
a + = 0.5/im, D ds = lOpc, a_ = 0.005/im, q = 4.5, and s = 3.0. 
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Fig. 4. — Spectral evolution of dust-scattered X-ray echoes. The solid line shows the initial 
Band spectrum with 5 = and E p = 200 keV. The spectra of the echoes at different times 
are shown from top to bottom. The dust parameters are a_ = 0.025/xm, a + = 0.25/xm, q = 
3.5, D ds = lOOpc, and s = 2.0. 




Fig. 5.— X-ray afterglows of GRB 060813 and GRB 060814. 5 w and E p w 200 keV are 
assumed. The consequent parameters are a + ~ 0.5 /xm and D ds ps 10 pc for GRB 060813 
and a + ~ 0.25 /im and -Dd s ~ 30 pc for GRB 060814. The other parameters, a_ ps 0.025 /zm, 
g ^ 4, and s ~ 2, are the same for both of them, [^ee the electronic edition of the Journal 
for a color version of this figure] 



